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We propose a mechanism of unconventional superconductivity in two-dimensional strongly- 
correlated electron systems. We consider a two-dimensional Kondo lattice system or double-exchange 
system with spin-orbit coupling arising from buckling of the plane. We show that a Chern-Simons 
term is induced for a gauge field describing the phase fluctuations of the localized spins. Through 
the induced Chern-Simons term, carriers behave like skyrmion excitations that lead to a destruc- 
tion mechanism of magnetic long-range order by carrier doping. After magnetic long-range order 
is destroyed by carrier doping, the Chern-Simons term plays a dominant role and the attractive 
interaction between skyrmions leads to unconventional superconductivity. For the case of the fer- 
romagnetic interaction between the localized spins, the symmetry of the Cooper pair is p-wave 
(p x ± ip y ). For the case of the antiferromagnetic interaction between the localized spins, the sym- 
metry of the Cooper pair is d-wave (d x 2_ y 2). Applications to various systems are discussed, in 
particular to the high-T c cuprates. 



I. INTRODUCTION 

Since the discovery of high-T c superconductivity in 
cuprates, jl| a large number of studies have been invested 
to uncover its mechanism of superconductivity. Although 
the issue is still in controversy, there are some specific 
properties concerning the mechanism of superconductiv- 
ity. First, the high-Tc cuprates show two-dimensional 
highly anisotropic behaviors. From their structure, the 
high-T c cuprates have a layered structure of CuO<i planes 
with several CuOi layers sandwiched between insulating 
layers. Reflecting this layered structure, measurements 
of the resistivity Q and optical conductivity ||J1] show 
strong anisotropic behaviors. Furthermore, the electro- 
magnetic properties of the superconducting state is well 
described by a Josephson-coupled layer model, m Sec- 
ondly, it seems that there is a close relationship between 
magnetism and superconductivity. In the absence of car- 
riers, the system is a charge-transfer insulator. Due to 
the large charge-transfer gap, a localized magnetic mo- 
ment is produced at each Cu site. These localized mag- 
netic moments form antiferromagnetic long-range order 
below Neel temperature. (Neel temperature is not equal 
to zero because of weak inter-layer coupling. ||) When 
holes are doped in the CuC>2 plane, they occupy O 2p a 
orbitals |9|-|l2|] and destroy antiferromagnetic long-range 
order. As we increase the hole concentration, Neel tem- 
perature decreases. Upon further doping, the antiferro- 
magnetic long-range order is destroyed and the system 
becomes the superconducting state. Also in the phase 
of superconductivity, antiferromagnetic correlations are 
observed. 

In contrast to conventional BCS superconductivity, 
superconductivity in the high-T c cuprates is unconven- 
tional. Symmetry of the Cooper pair is not s-wave but 
d-wave, or more precisely, d 2 .2_j / 2-wavc. ]l5| ] The fact that 
symmetry of the Cooper pair is d-wave, superconductiv- 
ity occurs in the vicinity of antiferromagnetic long-range 



order, and the absence of the isotope effect |16| suggests 
that the underlying mechanism of the high-T c cuprates 
be ascribed to the antiferromagnetic correlations. 

In addition, the structure of the C11O2 plane seems to 
play an important role for the pairing mechanism of high- 
Tc superconductivity. In La2- x Sr x CuOi system, sup- 
pression of superconductivity is observed at a structural 
phase transition point from an orthorhombic phase to 
a tetragonal phase. [ p7| Similar suppression of supercon- 
ductivity, which is induced by the same kind of structural 
phase transition, is also observed in La2- x Ba x CuOi sys- 
tem around x = 1/8. |L8| Since there is buckling of the 
C1LO2 plane in the orthorhombic phase, the fact that 
superconductivity occurs in the orthorhombic phase sug- 
gests that buckling of the C11O2 plane plays a significant 
role for the occurrence of superconductivity. 

In order to explain the mechanism of high-T c supercon- 
ductivity in the cuprates, a large number of theories have 
been proposed. Among others, the spin-fluctuation the- 
ory ]l9| proposes a pairing mechanism with d-wave sym- 
metry, or d x 2_ y 2 symmetry. For the d-p model, which is 
believed to capture the essential properties of the C11O2 
plane, the spin- fluctuation theory predicts d x 2_ y 2 pairing 
between d-orbital electrons at Cu sites. |^0|,^l| The same 
kind of d-wave pairing, i.e., d x 2_ y 2 wave pairing, is dis- 
cussed in a different context. In the resonating valence 
bond (RVB) theory [p2|p3| , spinons form d-wave pairing. 
pi| However, spinons are charge neutral quasiparticles 
and the electric current is carried by holons in the RVB 
theory. That is, d-wave pairing between spinons does not 
lead to d-wave superconductivity by itself. The d-wave 
pairing state of spinons rather describes the short-range 
antiferromagnetic correlations in the phase without anti- 
ferromagnetic long-range order. |25f| 

Although the relationship between these two d-wave 
pairings has not yet been clear, the fact that the doped 
holes occupy Oxygen p-orbitals ]9|-|Pz| and the sign of the 
Hall coefficient is positive pfj] indicates that carriers are 
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doped holes. We may say that it is doped holes to form 
the Cooper pair in the superconducting state. What we 
require to describe high-T c superconductivity is the d- 
wave pairing mechanism between doped holes. 

In this paper, we propose such a mechanism of super- 
conductivity. We consider a two-dimensional multi-band 
model which consists of a carrier system and a localized 
spin system with strong coupling between them. As a 
typical candidate of such a system, we consider a Kondo 
lattice system or a double-exchange system. The high- 
Tc cuprates can be described as a Kondo lattice system 
p7| , p8[ , where carriers are doped holes and localized spins 
are at Cu sites. In order to take into account the ef- 
fect of buckling of the plane, we consider spin-orbit cou- 
pling arising from buckling. We show that carriers in- 
duce frustration in the localized spin system in the pres- 
ence of spin-orbit coupling through Kondo or Hund cou- 
pling. This frustration effect can be described in terms 
of skyrmion excitations. The skyrmion excitation is cre- 
ated at each position of the carriers and plays a role of 
"magnetic" field for the carriers. Because of the "mag- 
netic" field produced around a carrier, the Lorentz force 
acts on another carrier. Due to this Lorentz force an at- 
tractive interaction is induced between carriers and leads 



to unconventional superconductivity. 29 3(J 

The outline of this paper is as follows. In Sec. [| we 
describe the model and the effect of spin-orbit coupling. 
In Sec. Ill , we show the mechanism of superconductivity 



in the case of ferromagnetic interaction between localized 
spins. Frustration effect induced by carriers is described 
as skyrmion excitations through a Chern-Simons term for 
the gauge field which describes the phase fluctuations of 
the localized spin system. The fact that carriers behave 
like skyrmion excitations provides a destruction mecha- 
nism of magnetic long-range order by carrier doping be- 
cause the magnetic long-range order is destroyed by the 
skyrmion excitations. After magnetic long-range order is 
destroyed by the skyrmion excitations, the Chern-Simons 
term plays a dominant role and the attractive interaction 
between skyrmions leads to p-wave superconductivity. In 
Sec. IV, we show the mechanism of superconductivity in 



the case of antiferromagnetic interaction between local- 
ized spins. In this case, the symmetry of the Cooper pair 
is d-wave, or more precisely d x 2_ y 2 wave. We also show 
that the doping carrier induces metal-insulator transition 
at the magnetic transition point. In Sec. |v|, we discuss 
applications to high-Tc superconductivity and other sys- 
tems. 



II. MODEL 



We consider a two-dimensional multi-band model 
which can be reduced to a model consisting of carriers 
and localized spins with strong coupling between them. 
Examples of such a model are the Kondo lattice system 
and the double-exchange system. In order to include the 



buckling effect, we introduce spin-orbit coupling arising 
from it. 



A. Hamiltonian 

The Hamiltonian of the Kondo-lattice system or the 
double-exchange system with the spin-orbit coupling 
term may be written in the following form: 

H = Hq + H int + H so + H spin . (2.1) 

Here the first term is the kinetic energy term for the car- 
rier system, 



H = -tY,(ch+H, 



(2.2) 



(id) 



where the summation is taken over the nearest neigh- 
bor sites. Carrier operators are represented by a spinor 
notation, 



( c h c ll 
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C 3l 



(2.3) 



The second term H mt represents Kondo or Hund cou- 
pling between the carrier spin Sj and the localized spin 



where Sj = ^c^acj with the components of er = 
(cti, (72,(73) being the Pauli spin matrices. We take 1, 
2, and 3 for the axes in spin space. We assum e tha t \J C \ 
is the largest energy scale in the Hamiltonian fl2.1| ). 

The third term H so in Eq. (2.1) represents spin-orbit 
coupling arising from buckling, 



(2.4) 



H* n = i 



E 



A^i "3 

3 77=(a,0),(0,o) 



crc 



j+V 



H.c. 



(2.5) 



where a is the lattice constant and the vectors 
(A^ , \^p\o) are proportional to both spin-orbit cou- 
pling of ions and the angle of buckling. The simplest 
example for the spin-orbit coupling term is presented in 
the appendix of Ref. pl[ | in which the spin-orbit cou- 
pling term for the s and p orbitals is derived. The spin- 
orbit couplingtcrms of the high-Tc cuprates are shown 



in Refs. [p2-34|. For simplicity, we assume the vec- 



tors for the orthorhombic phase of Y BaiCu-sOi-t,, 
that is, A(- a <°) - A( Q '°) and A^ = A<°'°> with 

| A (a,0)| = | A (0,a)| _ y 

The last term H sp i n in Eq. (|2.l[) represents the inter- 
action between the localized spins. For H sp i n we take the 
Heisenberg Hamiltonian 



H S pin — J ^ ' Si • Sj . 

(hi) 



(2.6) 
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In general, the Dzyaloshinskii-Moriya interaction is in- 
duced between the localized spins when there is spin- 
orbit coupling like H so . However, the Dzyaloshinskii- 
Moriya interaction does not play an important role in 
our mechanism of superconductivity. We neglect the 
Dzyaloshinskii-Moriya interaction in the following anal- 
ysis. 



B. Effect of spin-orbit coupling 

In the last subsection, we have introduced spin-orbit 
coupling arising from buckling of the plane. Since this 
spin-orbit coupling plays an essential role in our mech- 
anism of superconductivity, first we need to discuss the 
effect of spin-orbit coupling. 

The effect of the spin-orbit coupling term H so becomes 
apparent when we combine H so with the kinetic energy 
term Hq: 

Ho + h so = [ c t (- ta ° + lXm ■ °) c ^-n + Hc - 



i rj 



v^EE 

i T] 

(c. l+T1 + H.c] , 



(2.7) 



up to 0((\/t) 2 ) in the exponential, where <7o is the unit 
matrix in spin space. The factor exp (— | A^ ■ crj has 
the form of a unitary transformation of rotation in spin 
space. The axis of rotation is parallel to the vector A^ 
and the angle of rotation is 2\/t. 

The presence of this rotation at every hopping process 
of the carriers implies that the carriers introduce disor- 
der in the localized spin system through strong coupling, 
H lnt between the carriers and the localized spins. Disor- 
der produced by the carrier hopping processes provides 
a destruction mechanism of magneti c long -ran ge or der 



in the localized spin system. In Sees. [II D| and 



IV B we 



will show that this destruction mechanism of magnetic 
long-range order is represented as the effect of skyrmion 
excitations. 



III. FERROMAGNETIC CASE 

In order to illustrate the mechanism of superconductiv- 
ity, we first consider the case of ferromagnetic interaction 
between the localized spins, that is, J < in H spm . 



A. Schwinger bosons 

We are interested in the mechanism of superconduc- 
tivity based on the fluctuations of the localized spins. 
In order to describe the localized spins, we introduce 



Schwinger bosons. Description of the localized spin sys- 
tem in terms of the Schwinger bosons has some advan- 
tages. First of all, it is straightforward to describe the 
magnetic long-range ordered state. The phase with mag- 
netic long-range order is described by Bose-Einstcin con- 
densation of Schwinger bosons. J35||56| Another advan- 
tage is that we can directly construct rotation matrices 
for carrier's spins. Such matrices turn out to be useful 
for the description of the localized spin fluctuation effect 
on the carriers. 

Each localized spin can be described by Schwinger 
bosons, 



6 ii 



Z 3l 



(3.1) 



Here and Zj a are Schwinger bosons at site j and 



3° 



obey boson commutation relations: 



ij "aa' 



and 



n<7 ' ^jct' 



and the constraint 



J2 a z ja z ja = 25. ]35| , |3l| In the following, we consider the 
case of 5 = 1/2 for simplicity. However, it is straightfor- 
ward to extend the following arguments to general values 
of 5. 

In terms of Schwinger boson fields, the Hamltonian for 
the localized spin system can be written, up to constant 
term, in the following form: 



H s , 



where = 



\j\J2 flf - 

(ij) 



t 

ij * V ' 



(3.2) 



yt 



Turning to the path-integral for- 
malism, we introduce a Stratonovich-Hubbard field Q^ 
and to decouple the interaction term Fj.Fi 



ij ± n 



VzVzVX^VQVQ exp (-5 spin ) , (3.3) 



where 



5spin — 



cIt 



$>,(3 T -zAf) 

jo 



•J" 



(ij) L 



QnQ 



(Q 



where the t dependence of all fields is implicit and A^ B 
is introduced to impose the constraint z 



= 1. 



Now let us study the localized spin fluctuations. The 
spin fluctuations are represented by Q^ because we ob- 
tain (Qij) = *52 a (zi<jZj a ) at the saddle point. The spin 
fluctuation Qij consists of the phase fluctuations and the 
amplitude fluctuations. Since the latter turns out to be 
a high-energy mode, we focus on the phase fluctuations 
of Qij. 
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The phase fluctuations of Qij are connected with the 
local gauge transformation of ~Zj a (or z^ a ) and Zj a at 



U, 



each site. In fact, Eq. (3.2) is invariant under the local where 
gauge transformation Zj a — > exp(— i^ ). In the action 
S a -pim this gauge transformation involves a transforma- 
tion in the phase of Qij. That is, the transformation in 
the phase of can be described by a gauge field. Intro- 
ducing a gauge field and the mean amplitude Q = (Qij) 
and taking a continuum limit, we may write the action 
5 sp i n in the following form: 

5 spin = fdr frrV [z a (v, r) (d T - iA s T B - i\ SB ) z a (r, r) 
Jo J 

|J|Q I^ ^ „ ^ a l (3 5) 



cj — > UjCj, 




Cj Uj , 




( Z 3\ -Zjl \ 






z Jl \ 


V z 3i z it / ' 






z n J 



(3.9) 



(3.10) 



Under these transformations, H lni is reduced to H mt 
Hint — — x 

ing phase factor is introduced: 



4 c j a z c j- ■"■ n the hopping term, the follow- 



er j+rjUj = exp (—ir7 • Ajrj) 



(3.11) 



where 



^ B (r,r) = -i^^(r,r)<V CT (r,T). 



(3.6) 



Note that Eq. ( p.5[ ) is invariant under the gauge trans- 
formation: z a (r, t) — > £ CT (r, r) exp (— i#(r, r)), and 
^(r,r)^f (r,r)-3 M 0(r,T). 



B. Gauge field description of the strong correlations 

The spin fluctuations of the localized spin system affect 
the carrier system through strong coupling Hi nt - We may 
say that this strong coupling between the carriers and the 
localized spins is the origin of strong correlations. In or- 
der to take into account this strong correlation effect we 
rotate the spin of the carrier so as to be in the direc- 
tion of the localized spin at the same site. Through this 
transformation, the effect of the spin fluctuations on the 
carrier system is described by coupling to a gauge field. 

The action of the carrier system with 7f; n t is 



S c + S'u 



where 



dr 



^c,(r)(9 r - M ) Cj (r) + ff c + i/ II 



j V 



"j+r) 



(r)e- 



(3.7) 



^Ojir) 



-J c 5>(t).S,(t). 



(3.8) 



In order to rotate the carrier's spin Sj in the direction of 
the localized spin Sj, we perform the following unitary 
transformation: 



If the phase fluctuation Ajr/ is sufficiently slowly vary- 
ing, we can take the continuum limit. ( We w ill discuss 
the validity of this approximation in Sec. V A.) Thus we 
obtain 



S c + Si„t = dr d 2 r^(r, r^" 1 ^ + A^(r, r), 



(3.12) 

where = —id^ and A^ is the SU(2) gauge field arising 
from the spin fluctuations of the localized spins, 



An — 



\« / . - > 

ct=l,2,3 



A a u a a = -iUdJJ. 



(3.13) 



Note that A$, — *4?, B . The inverse of Green's function is 

H- 



G-\{K}) = (ik T + &) <r + g(k) • <t, (3.14) 

with £fc = k 2 /(2m) — \i (l/(2m) = t). Here g(k) is given 
by 



g(k) = 2\^k x + 2\ (0 ^k y - ^e 3 (3.15) 



As a result, the total action may be written in the fol- 
lowing form: 

S = S c + S ln t + >5 S pi n 

= f dT I d2r {^( r > T ) G ~ 1 ({^+-M)^(r,T) 
+ J2 M r * r) (d T - iAl - z a (r, r) 

a 

, \J\Q 



-\(V-iA 3 )z a (r,T) 



(3.16) 
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C. Effective action of the gauge field 



D. Skyrmion excitations 



From Eq. ( 3.16 ) we can see that the fluctuations of 
the localized spins affect the carrier system through the 
gauge field A^. Therefore, in order to investigate the 
effect of the spin fluctuations on the carrier system, we 
need to study the properties of the gauge field A^, that 
is, we need to calculate the effective action of the gauge 
field Aft- The effective action S_a, consists of two parts: 



Sa = S 



A 



Qspm 
A ■ 



(3.17) 



Here is the contribution from the carrier system and 
S s j£ m is that from the localized spin system. 

We obtain by integrating out the carrier fields. 
From Eqs. (3.14) and (3.15), one can see that A (r?) play 
a role of the Dirac 7 matrices in 2 + 1 dimension and J c 
plays a role of the Dirac fermion mass. The derivation 
of the effective action of the gauge field A^ is similar to 
that for massive Dirac fermions in 2+1 dimension. 37 4(J 
We find that the Chern-Simons term for the gauge field 
Afj, is induced. (Detail of the calculation is presented in 
the Appendix.) The induced Chern-Simons term has the 
following form: 



^4 ~~ 

Z7T 



r 

dr / d 2 vA\ (d x A 



(3.18) 



Here we retain only the third component of the SU (2) 
gauge field because it describes the ferromagnetic spin 
fluctuations of the localized spins. The coefficient of the 
Chern-Simons term is 



1 



9 = - sgn(J c A) 



with A = \[ aS)) \f a) - X { 2 a ' 0} \f' a) . Equation (|3J9|) is the 
expr ession at zero temperature. For finite temperature, 
Eq. (3.19) is slightly modified. However, we can neglect 
finite temperat ure e ffect as long as f3\J c \ ^ 1. 

The action (3.18) represents the combined effect of 
the spin-orbit coupling term H so and the strong cou- 
plin g ter m Hi nt . This effect is qualitatively described in 
Sec. II B. That is, destruction of magnetic long-range or- 
der. In the gauge field description, the effect is described 
by the Chern-Simons term for the gauge field. 

For the contribution from the localized spin system 
S*% in , it depends on whether there is magnetic long- 
range order or not. In the absence of magnetic long- 
range order, S s J^ ln may have a form of the Maxwell term. 
Since there is an extra derivative in the Maxwell term 
compared with the Chern-Simons term, we expect that 
the Maxwell term has unimportant effect for the long- 
wavelength and low-energy theory. On the other hand, in 
the presence of magnetic long-range order, the gauge field 
Af l (= A^) becomes massive since Schwinger bosons 
form Bosc-Einstein condensate. (That is, Schwinger 
bosons are in the Meissner phase with respect to the 
gauge field A^.) 



(3.19) 



Since coupling between the carriers and the localized 
spins is made only through the gauge field A^, we may 
write the effective action of the carrier system in the fol- 
lowing form: 



s: 



off 







dr / d ri/j(r,T) 



d T - fi + iA% 



i/)(t,t) 



+ £ J dT J drA * ( dxA y - 9yA ^ + sTn - (3 - 20) 

From this action, one can derive an important relation- 
ship between the carriers and excitations in the localized 
spin system. The variation of S*° ff with respect to A% 
yields 



rP<r(r,r) 



2?r 



ipxAy dyAj}) , 



(3.21) 



where sj = +1 and sj_ = — 1. This equation implies that 
a gauge flux is produced at the position of the carrier. 
This gauge flux corresponds to the skyrmion excitation 
similar to a topological excitation of the non-linear 
sigma model or CP 1 model p2]447| 

The gauge fluxes produced by each carrier play a role 
of vortices introduced in a BCS superconductor. In a 
BCS superconductor, which is Bosc-Einstein condensa- 
tion of Cooper pairs, disorder is introduced by vortices, 
or the electromagnetic gauge fluxes. As mentioned in 
Sec. [II A, the magnetic long-range order in the local- 



ized spin system is described by Bose-Einstein condensa- 
tion of Schwinger bosons. Since the Schwinger bosons 
couple to the gauge field At, its gauge fluxes intro- 
duce disorder in Bosc-Einstein condensation of Schwinger 
bosons. This implies that disorder is introduced by 
skyrmion excitations, or carrier doping through the rela- 
tionship (3.21). In the presence of magnetic long-range 



order, the skyrmion excitations have an excitation en- 
ergy gap. After magnetic long-range order is destroyed 
by skyrmion excitations, they become gapless excitations 
and the Chern-Simons term plays a dominant role in 
long-wavelength and low-energy physics. 

Before moving on to study the interaction between 
these skyrmion excitations, we make some comments on 
the time-reversal symmetry. From the coefficient of the 
Chern-Simons term Eq. (3.19), we see that the Berry 
phase induced by the Chern-Simons gauge flux is 2tt. 
Such a Berry phase preserves the statistics of particles. 
By contrast, in the anyon system |^,|^| the Berry phase 
is 7r/2. The Berry phase of non-integer multiples of tt im- 
plies the time-reversal symmetry breaking, ]50|| whereas 
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in our case, there is no such implication of the time- 
reversal symmetry breaking arising from the Berry phase. 
In addition, there is no mean "magnetic" field as long 
as s a p a {r, t)) = 0. On the other hand, since the 
gauge field describes the localized spin fluctuations, both 
the time-reversal symmetry and the parity are broken in 
the localized spin system by the presence of the Chern- 
Simons term. 



E. Superconductivity 

In the phase without magnetic long-range order in the 
localized spin system, the Chern-Simons term plays a 
dominant role in the action Sa- I n this phase, an attrac- 
tive interaction is induced between skyrmions. Through 
this attractive interaction, carriers form the Cooper pair. 

Intuitively, we can understand the occurrence of an 
attractive interaction between skyrmions as follows. As 
discussed in Sec. pi B| , the carrier rotates its spin at each 
hopping process due to the spin orbit coupling term _ff so . 
These rotations of the carrier's spin affects the localized 
spins through the strong coupling term H mt . This effect 
can be described as the formation of a spin configuration 
in the localized spin system. This spin configuration car- 
ried by each carrier generates a "magnetic" field around 
other carriers through the Berry phase. Therefore when 
a carrier passes another carrier with velocity, a Lorentz 
force acts between them. This Lorentz force plays a role 
of pairing interaction. 

Now let us go into detail. The Lorentz force is de- 
rived from minimal coupling between the carriers and 



the gauge field A 



Mnt = / d 2 r^- ]T s^l{v) (k • A 3 + A 3 ■ k) V^v(r). 

£7 

(3.22) 



From Eqs. ( |3.2l| ) and fl3.22| ), we eliminate the gauge field 
A 3 upon taking the Coulomb gauge. Thus we obtain 



1 

2Q 



H = J2&ci a c ka + ^ s ai s a2 v klk2 

ka" ki^k 2 crier 2 

XC kieri C -kiCT2 C - k 2 cr 2 C k2CTl • (3.25) 

Here we set q = in the interaction term to focus on 
the equilibrium state. We investigate the po ssibil ity of 
superconductivity based on the Hamiltonian ( 3.25 ). For 
simplicity, we analyze the Hamiltonian within a mean 
field theory. We introduce the following pairing matri- 
ces: 



xc 



-ki+q/2,CTi L --ki+q/2,cr2 C - k 2+q/2,CT2 C k 2 +q/2,(T 1 

where fi is the area of the system and 
4tvi ki x k 2 



(3.23) 



m6 Iki — k a l 



(3.24) 



Equation ( 3.23 ) represents the interaction between the 
carriers mediated by the gauge field A 3 ^. Note that in 
this equation no parameters characterizing the skyrmion 
excitations appears except for the Chern-Simons term 
coefficient 6. Furthermore, there is no retardation effect. 

Now we may write the effective Hamiltonian for the 
carrier system in the following form: 



A ^2 = ^ £ Vkk'^-k'^Ck'^), 



(3.26) 



k'(^k) 



k'(#k) 

In terms of these matrices, we define the mean field 
Hamiltonian 



krr 



ktTiO-2 



H M F = (& C L C k<T - &c- kCT c T _ kCT 

T 

£ ' [( Ak )L c -- c - + ^AJ-^ 

(3.28) 



Here the summation in k space is taken over a half of the 
first Brillouin zone. 

For the spin singlet pairing case in which A k | = 
— A k T = Ak, the gap equation can be derived by tak- 
ing the standard procedure: J5lj 

1 v-^ A k , 8E k , 

k'(#k) 



A, 



with £ k = v^fc + [A 



with A k T 



A k 



and A k ^ 



For the spin triplet pairing case 
we obtain the same 



A k 



gap equation ( 3.2E ). A pairing state with A^ ^ and/or 



A ii ^ ^ ma y be stabilized in the presence of an in-plane 
magnetic field or at the sample's boundary. Here we do 
not consider such a c ase. 

The gap equation ( |3.2S| ) is the same as the gap equa- 
tion for the composite fermion pairing state at half-filled 
Landau levels. |5^,|53| We apply the same analysis of 
Rcf. [ [52] ]. We assume that the gap function Ak has the 
following form: B3] 



A k = A fc exp (-ii&k) 



(3.30) 



After substituting Eq. ( |3.30|) into Eq. ( ^.29| ), we integrate 
the angular variable #k using Cauchy's theorem by tak- 
ing exp(i£$k) as a complex variable. In this calculation, 
we find that the attractive interaction arises only in the 
case of I > 0. From the fact that the case of £ < yields 



G 



a repulsive interaction, we may exclude the possibility of 
a linear com bination of components with i > and £ < 
in Eq. ( pC| ). 

For the ground state, we obtain |52l 



A fe = 



2m9 



dk> k -^ 
k Ek> 



k 
A 7 

(3.31) 



In order to solve this non-linear integral equation, we use 
an approximation. Fro m the asymptotic behavior of the 
right hand side of Eq. fl3.3l| ), we set 



f AeF(k/kpY, for k < k F} 
\ AeF(kF/k) e , for k > kF- 



(3.32) 



Substituting this equation into Eq. (3.31), we obtain |5 



{dk/ki 



(k/k F ) 



21+1 



[(k/k F ) 



(dk/kp) 

■l[{k/k F ) 2 



-1] + 
(k/k F 



A 2 (k/k F y e 

1-21 
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A 2 (k/k F )~ 2i 



The remaining parameter A can be evaluated numeri- 
cally from Eq. (3.33) . |p3|l For each £ we estimate the 
left hand side of Eq. ( 3.33 | ) w ith varying the value of A. 
The point at which Eq. ( 3.33j ) is satisfied gives the value 
of A. From this analysis we find that the largest gap 
is obtained for the case of £ = 1 and A e=1 ~ 3.1. Fur- 
thermore, this state has the lowest ground state enegy. 
Therefore, the grou nd state is p-wave superconductiv- 
ity. From Eq. ( |3.30| ), the symmetry of the Cooper pair 
is p x ± ip y . Incidentally, this is the same pairing state 
as that proposed in Sr 2 RuOn. |54| ] However, we cannot 
apply our pairing mechanism to this system. We shall 
discuss this point in Sec. |v| 



IV. ANTIFERROMAGNETIC CASE 



c, — > UiCi (i S A) 
Cj -> Uj (-ia 2 ) Cj (j e B) 



(4.1) 



The definition of Ug is given by Eq. (3.10). Note that the 
matrix, 



Uj (-ia 2 ) 



-_ z n -*jt 
z j1 -Zjl 



(4.2) 



can be derived from the matrix Uj by the transformation 



-jn 

Z 3l- 



(4.3) 



If we do the same transformation in S 7 - = iz 



tain S., 



a 1 



crz. Thus, in Eq. (|4.1| 



2^j<tz, we ob- 
the presence of 



the factor (— ia 2 ) at the B sublattice implies that the car- 
rier's spin is in the direction of the staggered component 
of the localized spins. 

In order to eliminate the factor (— ia 2 ), we perform the 
following transformation at the B sublattice: 



6. (3.33) 



(i e B). 



By this transformation, we obtain 



(4.4) 



S c + S int = I dr / d 2 ri;(r,T)G' 1 ({k^ + A f ,})^(r,T), 

(4.5) 



where G 1 {{k fl }) is given by Eq. (3.14). This action has 
the same form as the action of the ferromagnetic case, 
Eq. Q3.12| ). Therefore, the same Chern-Simons term is in- 
duced by integrating out the carrier fields. However, we 



need to perform the inverse transformation of Eq. (4.4) 
when we study the symmetry of the Cooper pair because 
the transformation (4.4) affects the order parameter of 
the Cooper pair. Furthermore, the action of the local- 
ized spin system, of course, differs from the action of the 
ferromagnetic case. In particular, the relevant gauge field 
component is different from that case. 



Now we study the mechanism of the case of antiferro- 
magnetic coupling between the localized spins. Although 
we can apply a similar analysis, the antifcrromagnetic 
case is more complicated than the ferromagnetic case 
because we need to cope with staggered magnetization. 
That is, we need to distinguish the A and B sublattices. 
Fortunately, there is a transformation by which the sys- 
tem is mapped onto a similar model of Eq. ( 3.16 ) in 
Sec. HI. We introduce such a transformation and an- 



alyze the mechanism of superconductivity through the 
transformation. 

As in the ferromagnetic case, we introduce Schwinger 
bosons and rotate the carrier's spins so that each of them 
is in the direction of the localized spin at the same site. 
The rotations are performed by the following transfor- 
mations: 



A. Action of the localized spin system 

In order to identify which component of is con- 
nected with the antifcrromagnetic spin fluctuations, we 
calculate the action of the localized spin system. In terms 
of the Schwinger bosons, the Hamiltonian is written as, 
up to constant, 

-ffspin = 77 / ] z i<j 1 z i^2 Z i+rj,a2 Zi + r l- a 'i " 



ieA T]=(±a,0),(0,±a) ci,<t 2 



(4.6) 



To write this Hamiltonian in a tractable way, we perform 
the following transformation at the B sublattice p3] 
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Z 3t 



(4.7) 



Turning to the path-integral formalism, we obtain 
-Z S pin = J T>~zDzT>\ SB exp(-S , spin ), where 



dr 



EE^+ A f) 



■EEE" 

ieA 7] <T 1 (T2 



\-Tf,ai Z>i<7\ ^i-\-Tf,a2 



(4.8) 



where the r dependence of all fields is implicit. In 
order to decouple the interaction term, we introduce 
Stratonovich-Hubbard fields Qi.i+r) and Q, L 2 sp i n = 
/ VzVzVX SB VQVQ exp (-5 spin ), 



dr 



EE^- (d T + \ s e B )zea 



E E f - < < • r/ ' - '/ 



-iEEE(^ 



A+rjZi+ri.a* 



(4.9) 



The spin fluctuation field Qi^+ri consists of the phase 
fluctuations and the amplitude fluctuations. The latter 
is irrelevant for our analysis as it is in Sec. III. The 
phase fluctuations are connected with a gauge invariancc 
of Schwinger bosons. |35| We include these phase fluc- 
tuation degrees of freedom later by imposing the gauge 
invariance. 

We set Qi ti+ r] = \Qi,i+r]\ = Q = const, and Af B = 
Asb = const.. Then, the action is diagonalized in k- 
space. Introducing the following fields, 

Ckff = ^ [(Zkcr + Zk+Q,o-) + (Z-kcr - Z-k+Q,a)} , (4.10) 



5 kc r = - [(z ka + Zk+Q, CT ) - (Z-ka - 2-k+Q,<j)] , (4.11) 

with Q = (jr/a,Tr/a), the action is written in terms of 
these fields as 

Sspin = E E E [~ iUn (^ko-Cko- + Ckcr^ko-) 

iuj n k u 



the mass of £ ko - is identically zero in the ordered phase. 
Therefore, we can safely integrate out S k(T and we obtain 



Sspin — EE 

k a 



AsB — £k 



AsB + £k 



CtoCk*. (4.13) 



Taking the continuum limit and recovering the gauge 
invariance of the Schwinger bosons, we obtain 

Sspin = - r^dxo [ rf 2 xE [l (9, + Uf ) C<t| 2 + 

9 Jo J „ L C sw 



(4.14) 



where g = 4V2a, A sw = V A Ib = 4J 2 Q 2 , c sw = V2JQa, 
and Xq = c sw t. In these pa rame ters we set Asb = 2 JQ 
except for A sw . Equation (4.14) is invariant under the 
gauge transformation [pq] 



C(r, r) -> C(r, r) exp [i0(r, r)] , 



(4.15) 



^B(r,r)^^B(r,r)-9^(r,r). (4.16) 
This gauge transformation corresponds to 



Ziexp(i9i) (ieA), 
Zjexp(-i0j) (j e B), 



(4.17) 



because if we take the set of (even, even) and (odd, odd) 
for the A sublattice, then 



Qa = Zt a for IeA, 

Civ = zt<r for £ G B. 



This equation is verified as follows 

1 v-v 

= j\rE ^ kCT ex P( ik ' 

k 



(4.18) 



-E 

k 



z k(T e ikRf (l + c- 4Q Rf ) 



2N^ 

k 



l-e 



-iQR, 



(4.19) 



Since the gauge field is connected with the phase fluc- 
tuations of Qi y i+r), the gauge field ^ B has the following 
form 



A 



SB 



(4.20) 



In order to find the relationship between _4^ B and A^, 



+ (A SB + £k)C kCT Ck CT + (Asb - e k )S kff S kCT ] , (4.12) we ™} te A* ™ terms of Q a and Q G . Thus, we find A 



where e k = — 2 JQ[cos(k x a) + cos(k y a)] and the summa- 
tion in k space is taken over a half of the first Brillouin 
zone. One can see that the mass of S k(T is Asb + 2 JQ, 
which is larger than the mass of Cka-, Asb — 2JQ. |Hj| 
Furthermore, the mass of S k(T is non- vanishing whereas 



*4 SB . From the gauge invariance of the Schwinger 
bosons, one can see that there is a correspondence be- 
tween A\ and .4^ B . Therefore, A^ is connected with the 
antiferromagnetic spin fluctuations. 

As a result, we may write the effective action in the 
following form 



8 



S — S c + S sp i n + Scs 



C. Superconductivity 



dr J d 2 ri/i(r, r) 
d T + iAl - n + ^- (-*V + A 1 ) 2 



2^ 7o 



$(r, r) 



A? 



\(d»+iAl)Ux)\ 2 + -^\Ux)\ 



dr / d 2 r^ (d x A\ - d y Al) 



(4.21) 



Now we investigate the property of superconductivity. 
In order to identify the symmetry of the Cooper pair, we 
need to perform the inverse transformation of Eq. ([Oh. 
However, since the calculation of the pairing matrix is 
much easier for the system after the transformation ([Oh, 
we investigate the pairing matrix through the transfor- 
mation (4.4). If the B sublattice is consists of the sites 
with (even, odd) or (odd, even), then the transformation 
(4.4) is written in k-space as 



Ck - c k +Q = i<7 2 (2k - Ck+Cj) . 



(4.23) 



B. Skyrmion excitations 



Since there is no change at the A sublattice, that is, 
c k + Ck+Q = c k + c k+ Q , we obtain 



As discussed in Sec. HID, carriers are connected with 



skyrmion excitations in the localized spin system through 
the Chern-Simons term. However, the connection is 
slightly different from the ferromagnetic case because the 
spin fluctuations are described by the gauge field Az in- 
stead of A^. If we take the 1-axis for the quantization 
axis in spin space, the relationship between the carrier 
and the skyrmion excitation is 



s CT Vv(r,T)'!/v(r, r) 



2tt 



{d x A l y -d y Al). (4.22) 



Contrary to the ferromagnetic case, a significant fea- 
ture appears for the antiferromagnetic case, that is, a pin- 
ning mechanism of carriers in the antiferromagnetic long- 
range ordere d pha se. This can be seen as follows. The 
r elat ionship (4.22) is obtained after the transformation 
(4.4). In order to capture the proper nature of skyrmion 
excitations, we mus t go back to the frame before the 
transformation (|4.4| ). Performing the inverse transfor- 
mation of Eq. (|4.4[) at the B sublattice, we find that an 
additional si gn ch ange is brought about in the left hand 
side of Eq. ( 4.22 ), that is, a skyrmion (anti-skyrmion) 
excitation transformed into an antiskyrmion (skyrmion) 
excitation. Therefore skyrmions or antiskyrmions can- 
not move to the nearest neighbor sites as long as there is 
antiferromagnetic long-range order and skyrmion excita- 
tions have a gap. This suggests an insulating behavior of 
the carriers in the antiferromagnetic long-range ordered 
phase. 

Although this is a new feature which appears in the an- 
tiferromagnetic case, the destruction mechanism of mag- 
netic long-range order is the same as that discussed in 
Sec. [II D. Antiferromagnetic long-range order is de- 



stroyed by carrier doping because carriers behave like 
skyrmion excitations. After antiferromagnetic long-range 
order is destroyed, skyrmion excitations become gapless 
excitations. The Chern-Simons term plays a dominant 
role and the attractive interaction between skyrmions 
leads to superconductivity. 



c k = - [{(T + «(T 2 )c k + (er - ^ 2 )c k+Q ] , (4.24) 
Ck+Q = g [(°0 - «cT2)5k + {(To + ^Sic+q] . (4.25) 



In order to apply the analysis done in Sec. HIE 
need to introduce the following fields 



we 



c k | 



1 



V2 U -1 



(4.26) 



because in Eq. (4.22) we take the 1-axis for the quantiza- 
tion axis in spin space. In terms of these fields Xko-, the 
spin singlet order parameter is written as 



(c_kj.c kT - c_ kT c ki ) 

= 2 Kx-klXkT - X-kTXki) 

+ (X-k+QJ.Xk+QT - X-k+QTXk+QJ.! 



(4.27) 



For the spin singlet pairing state, the gap function A k 
is given by 



A k =fA( 1) +A^)/2, 



(4.28) 



where A^ 1 ' is the gap function for (x- k |Xkt) = 
-(X-kTXki) and A^' is that for (x- k +QiX k +QT) = 
— (X-k+QtXk+Qi)- In the continuum limit, A k satis- 
fies, at T = 0, 



A 



(i) 



1 A (1) 

2n E ykk '^T)' 
k'(#k) ^k' 



(4.29) 



where E 1 
satisfies 



(i) 



H + lA^I 2 . On the other hand, A 



(2) 



A k 2) = 



1 



X(2) 

E 



k'(#k) 



(2) ' 



(4.30) 
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where E,„ 



^ - \j^i + \^k\ 2 - ine mmus si s n 

in the right hand side of Eq. (4.30) comes from 



The 



the sign change in minimal coupling to the gauge 
field. The minimal coupling term is derived from 

G^EXik + M)U^o,aA, where G£ E X{k M }) 



represents the kinetic energy part of the inverse of 
Green's function G _1 ({fc M }). In the lattice system, the 
displacement of Q in k space yields the sign change 
because, — 2t [cos (k x + Q x ) a + cos (k y + Q y ) a] = 
+2t (cos k x a + cos k y a). This sign change affects the 
interaction term through the minimal coupling term. 
Equation ( 4.3C ) is understood to be the equation ob- 
tained after taking the c ontinu um li mit. 

The analy sis of Eqs. (|4.29|) and ([4.3 0| ) parallels that 
in Sec. HIE. So we also take Eq. (3.30) for the form of 



A k 1} and A k " ; . From Eqs. fl4.29p and (|4.3C|), we sec that 



(2) 



A k ^ and A^ J have opposite chirality. For the case of the 
spin singlet pairing, \£\ takes the val ues \£ \ = 2,4,6, 



From the same analysis done in Sec. HIE, we find that 
the pairing state with the lowest energy is \£\ = 2, or 
d-wave. The pairing function is given by 



A k = A fe [exp (2i0 k ) + exp (-2i0 k )] /2 
= Afe (cos 2 k - sin 2 6» k ) , 



(4.31) 



where Afe = Aep(k < /k > ) 2 with A ~ 1.3. Here if k is 
smaller than kF then k < = k and fc> = kp and vice 
versa. From Eq. ( 4.31 ), we see that the precise symme- 
try of the Cooper pair is d x i_ y i. 

For the spin triplet pairing case, we find that A^ = 
A = and A|| = Aj^ cx sin0 k . In the d- vector no- 
tation of the triplet pairing state, ]56] ] this is written as 
dk = k y e^. From the analysis of the Ginzburg-Landau 
free energy within the two-dimensional representa- 
tion in which the basis states are {k x es, k y ez}, one can 
see that the free energy of this pairing state is higher 
than that of d k = [k x ± ik y )e^. This suggests that the 
pairing state with d k = k y e^ is not stabilized and we can 
exclude the possiblity of the spin triplet pairing state. 

From the above analysis, we may conclude that the 
symmetry of the Cooper pair is d-wave, or more precisely 
d x 2_ y 2-wsLve. This is the same pairing state as that of the 
high-T c superconductors. 



V. DISCUSSION 

In this section, we discuss the condition for that the 
Chern-Simons term is induced and applications to the 
high-T c cuprates and other systems. 



A. Condition for the Chern-Simons term 

The derivation of the Chern-Simons term is based on 
the continuum approximation. This approximation is 



justified when the length scale of the gauge field is rel- 
atively larger than the length scale determined by spin- 
orbit coupling. We first discuss the validity of this as- 
sumption. 

As we have shown in Sec. 



II B 



the effect of spin-orbit 
coupling is to rotate the carrier's spin at every hopping 
process. The angle of this rotation is 2A/t(= 9 so ) for 
the nearest-neighbor hopping. On the other hand, the 
length scale of the gauge field is determined by t he fluc- 
tuations of the spin system through Eq. (3.13). This 
length scale is given by hc sw / (7rA sw ) from the analysis of 
the Schwinger bosons and it is translated into the fluctu- 
ation angle of the localized spins. For the nearest neigh- 
bor sites, the fluctuation angle is, ~ a/(hc sw /irA sw ) = 
it A sw / (hc sw I a) . In order to apply the continuum approx- 
imation to the gauge field, this fluctuation angle should 
be smaller than 9 so = 2X/t, that is, 



2A 



ttA s 



t ^ (hc BW /a) 



(5.1) 



This is the condition of taking the continuum limit for 
the gauge field. 

The condition ( |5.l| ) is satisfied as long as sufficient 
magnetic correlations are preserved, or the gap of spin 
wave e xcita tions is relatively small. Note that the con- 
dition (5.1) is just for the localized spin system. It only 
provides the condition for the existence of the Chern- 
Simons term. Although our mechanism of superconduc- 
tivity is based on the presence of the Chern-Simons term, 
the superconducting state does not rely on directly the 
value of spin-orbit coupling A. In fact, the gap of super- 
conductivity is independent of the parameter A as shown 
in Sees. Ell and FVCl. 



B. High-T c cuprates 

As we have discussed in Introduction, the high-T c 
cuprates can be described as the Kondo lattice system. 
In addition, buckling seems to play an important role 
concerning the occurrence of superconductivity. If buck- 
ling and the parameters characterizing the localized spin 
system fulfill the condition (5.1), then the Chern-Simons 
term is induced. 

If we can apply our mechanism to the high-T c cuprates, 
then it describes the underdoped region. Because the un- 
derdoped region is close to the antiferromagnetic long- 
range ordered phase. Such a phase is properly described 
by the Schwinger bosons. However, in the optimal doped 
region and the overdoped region of the high-T c cuprates 
the localized spin system is much more disordered than 
in the underdoped region. 

In order to describe strongly disordered regions, the 
Schwinger bosons are not appropriate fields. The de- 
scription in terms of fcrmionic fields is more suitable than 
the bosonic description of the localized spin system. The 
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description by fermionic fields instead of the Schwinger 
bosons corresponds to describe the Cu site degrees of 
freedom by fermion fields. Such fermionic degrees of free- 
dom may be observed by angle-resolved photoemission 
spectroscopy experiments (ARPES). 

It should be noted that the antiferromagnetic corre- 
lation described in terms of fermion fields has the form 
of the singlet pairing between the fermions, because the 
order parameter for the antiferromagnetic correlations, 
Qi,i+r) in the Schwinger boson system is replaced by 
(/»T fjl ~ fil fit) = fij i in terms of the fermion fields, f ia ■ 
In deriving this relationship we need to go bac k to the 
original system before the transformation of Eq. (4.7) and 
we use the constraint f} a fi a — 1- This correlation is 
similar to that characterizes the RVB state. |25| From the 
mean field analysis, we find that shows d x 2_ y 2 sym- 
metry. J24[] Although this d x 2_ y 2 symmetry is the same 
as that of the Cooper pair of holes, they have completely 
different origins and must be discussed independently. 

Now we discuss some properties of the superconduct- 
ing state. In our mechanism, there is a ch aract eristic 
excitation. Because of the transformation ( 4.23 ), two 
components of Cooper pairs appears: (x-kj,Xk|) an d 
(X-k+QiXk+Qt)- Therefore, there is an excitation be- 
tween them that creates one quasi-hole in the one com- 
ponent of the Cooper pairs and one quasiparticle in the 
other with the momentum Q = (tt/o, ir/a) and the en- 
ergy 2 A kF . Such an excitation can be detected by inelas- 
tic neutron scattering. This excitation may be identified 
with the 41meV peak observed by the neutron experi- 
ments. |37j 

In addition, it should be stressed that the strength of 
spin-orbit coupling does not affect the transition tem- 
perature of superconductivity as long as the condition 
for the presence of the Chcrn-Simons term is satisfied, as 
discussed in the last subsection. 

Buckling affects superconductivity not through the 
spin-orbit coupling term but through the Fermi energy 
because the superconducting gap is proportional to the 
Fermi energy. If we increase the angle of buckling, then 
the hopping amplitude t in the plane may be reduced. 
Such a reduction leads to decrease of the Fermi energy. 
Therefore, if we increase the angle of buckling, the tran- 
sition temperature is rather reduced. This is consistent 
with the experiments. |58j 



our mechanism. Secondly, we require the system with the 
layered structure and buckling of the planes. Although 
La,2-2xSri + 2 X Mn20r is a layered double-exchange sys- 
tem, the x = 0.3 compound shows ferromagnetic correla- 
tion. |5SJ If the compound with buckling and small x in 
which the sample shows antiferromagnetic correlation is 
provided, d x 2_ y 2 superconductivity based on our mecha- 
nism will be realized. 



D. Other systems 



From the symmetry of the Cooper pair in Sec. [II E , 
one might think of the application to SY2-RUO4. How- 
ever, the magnetism relevant to this system is itinerant 
ferromagnetism. In fact, all of the relevant d-orbitals 
d X y, dyz, and d zx form conduction bands and the Fermi 
surface of them are observed by the de Haas- van Alphcn 
effect. |6(| In such a system, we cannot expect the forma- 
tion of the localized spins. Therefore, we cannot apply 
our mechanism to S^RuOi. 



For the application of the mechanism of Sec. Ill to real 



materials, we require a ferromagnetic superexchange in- 
teraction between the localized spins. In order to produce 
such a superexchange interaction, we need at least three 
kinds of ions or multi-band structure for the magnetic 
ions to constitute the conduction layers. Furthermore, of 
course, we need buckling of the planes. 

At present, it seems that there is no material with all 
of these properties. However, if such a system exists, 
then we can expect higher superconductiving transition 
temperature than the high-T c cuprates. Because the gap 
is larger than that of the antiferromagnetic case within 
the analysis in which the long-range Coulomb interaction 
is neglected. 
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APPENDIX: INDUCED CHERN-SIMONS TERM 



C. Double-exchange systems 

Our mechanism can be applied to the double-exchange 
system. However, in application to that system the fol- 
lowing conditions should be satisfied. First, we must 
detect superconductivity in the region where the an- 
tiferromagnetic correlation between the core spins are 
preserved. That is, the carrier number must be small. 
Whereas the region in which the ferromagnetic correla- 
tion between the core spins dominates, we cannot apply 



In this appendix, we derive the effective action for the 
gauge field arising from the carrier sys tem. Integrating 
out the carrier fields from the action ( 3.12 ), we obtain 
the effective action S' 



A- 



S C A = -Tr \ In 



{in [G- 1 ^ + A,)] - ln [G- 1 ^)] } ■ ( A1 ) 



We expand this action with respect to as S A 
+ S (3) + •••, where 
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S< 2 > = -Tr 



5 (3) 



-Tr 



1 /9G- 1 (dG- 1 A 



(A2) 



G 2 \ dk x ' Av 



In order to calculate S^ 2 \ we apply the derivative 
expansion technique. To illustrate this technique, let 
us consider a simple one dimensional example: J = 

Tr Vi{k)F 1 (x)V 2 (k)F 2 {x) , where k = -id x and Vj (k) 

and Fj(x) are functions of k and x, respectively. Insert- 
ing the identities J dx\x)(x\ = 1 and J2 k \fy(k\ = 1) we 
obtain 

J= I dxJ2(x\Vi(k)F 1 {x)\k)V 2 (k)F 2 {x)(k\x) 
J k 

= / dxY^F^v^y^^-p-ix^F^xmik^) 

dx^F 2 {x)V 2 {k) 

k 

v ^P- v "! ( ~* r ,M m \x)k n - m (x\k)(k\ X ) 

^ ra ^ m!(n-m)! 1 

n=0 m=0 v y 



I^V^V^fe) j dxF 1 {x)F 2 {, 



k 

-Y 



V{(k)V 2 {k) / dxF[{x)F 2 {x) 



+Vx{k)V 2 \k) J dxFi(x)F%(x) 
-^tJ2 V "W V *W [ dxF{'(x)F 2 (x) 



(A4) 



In the second term in the last line we have taken into ac- 
count the term obtained by the partial integral. Applying 
this derivative expansion and retaining the most domi- 
nant term in long- wavelength and low-energy physics, the 
action S^ 2 ' is evaluated as 



^ (2) = i E E I dr I d 2 vA«d„A a x , (A5) 



where 

jaa = n 



tr 



<9fc„ 



(A6) 



Here the trace is taken over spin space. All functions 
I/luX can be calculated in a similar way. For example, 



(*<■■ 



0) A (0,a) _ A (a,0) A (0,a)^ ^ 



(A7) 



In deriving this equation, we have taken into account 

(77) 

another choice of X a ' connecting coordinate space with 



spin space. In Eq. (A7), K is given by 



(A3) 



n 4-^ 



i 



, Pv+&) 2 -|g(k)i 2 r 
rf 2 k r/'(a + |g(k)i)+/'(a-|g(k)i) 



(2tt) S 



4|g(k)] s 



/(& - |g(k)|) - /(& + |g(k)|) 



4|g(k)|3 

In the limit of /3|J C | — > co, this equation is reduced to 



(A8) 



K = 



1 



87rA|J r 



■ tanh 



P\Jc 



(A9) 



Thus, S^ 2 ' has the form of the Chern-Simons term, 

+A 2 T (d x A 2 y - d y Al)] . (A10) 
The coefficient of the Chern-Simons term is given by 

0=isgn(J c A)tanhS j (All) 
A o 

withA = Ah 0) A^ Q) -A^ 0) A^ Q) . 

From similar calculations, we find that has a form 
of a non-Abelian Chern-Simons term, 

S ( cl M r ( A * A l ~ A * A l) ■ ( A12 ) 

This non-Abelian Chern-Simons term can be reduced to 
an Abelian Chern-Simons term upon using the curl free 
condition, 2{A\A\ - A 2 x A y ) = d x A\ - d y A%: 



? (3) 

7 CS 



2^ Jo 



dr / d 2 vA\ (d x A 3 y - d y Al) (A13) 



As a result, we may write the effective action for the 
gauge field arising form the carrier system in the follow- 
ing form: 



-A\{d x A\ - d y Al) - A 2 T (d x A 2 y - d v A 2 x )] ■ (A14) 
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